Existence of Dynamical Scaling in the Temporal Signal of Time Projection Chamber 
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The temporal signals from a large gas detector may show dynamical scaling due to many correlated 
space points created by the charged particles while passing through the tracking medium. This has 
been demonstrated through simulation using realistic parameters of a Time Projection Chamber 
(TPC) being fabricated to be used in ALICE collider experiment at CERN. An interesting aspect 
of this dynamical behavior is the existence of an universal scaling which does not depend on the 
multiplicity of the collision. This aspect can be utilised further to study physics at the device level 
and also for the online monitoring of certain physical observables including electronics noise which 
are a few crucial parameters for the optimal TPC performance. 



Dynamical scaling refers here to powerlaw distribution 
or a powerlaw with an exponential cut-off that describes 
certain correlation phenomena in the dynamical systems 
distinctly being different from the random statistical pro- 
cesses. In the thermodynamical context, it describes a 
critical phenomena associated with the phase transition 
[1] while in many other complex systems it corresponds 
to the so called self organised criticality [2]. Although a 



powerlaw distribution (/ 



) is a common dynam- 



ical feature, a powerlaw with an additional exponential 
(/ ~ x~^e~ ax ) makes the distribution normalizable for 
all values of (3 (where a and (3 are positive constants) and 
also many real world systems like World Wide Web and 
social networks show this cut-off [3]. In this paper, we 
show the presence of both type of scalings, a powerlaw 
and a powerlaw with exponential in the temporal signals 
(comprising of time gap and bunch length distributions) 
of a large Time Projection Chamber (TPC) which is a 
type of gas detector used for three dimensional track- 
ing of the charged particles passing through it during a 
particle physics experiment. The time gap distribution 
along the drift direction shows a dynamical scaling which 
is independent of the multiplicity of the collisions above 
a critical value. While this scaling behavior itself is an 
interesting aspect to investigate the correlation phenom- 
ena in a gas detector at the device level, it can be further 
utilised for online monitoring of certain important TPC 
parameters including the constant rise in electronics noise 
level due to long exposure to radiation. 

The TPC is the main tracking device of the ALICE (A 
Large Ion Collider Experiment) [4] at the Large Hadron 
Collider (LHC) at CERN optimized for the study of 
heavy ion collisions at a centre of mass energy ~ 5.5 
ATeV. It is a large gas filled detector of cylindrical de- 
sign with an inner radius of about 80 cm, an outer radius 
of about 250 cm, and an over all length in the beam direc- 
tion of 500 cm. A charged particle passing through the 
gas volume creates electrons by ionization. The electrons 
drift in the electric field towards the read out chambers 
(multiwire proportional counters with more than 550000 
cathode pads read out located at the two end-caps of the 
TPC cylinder) where they are amplified in the field of the 
sense wires by a factor of several 10 4 . This signal is cou- 
pled to the read out pads which are on ground potential 



and at a few milimeters distance behind the sense wire- 
plane. The detail aspect of the design and simulation re- 
sults can be found in ref [5] . For simulation, the charged 
particles (mostly pions) with different multiplicities are 
generated using HIJING parametrization corresponding 
to Pb + Pb collisons at 5.5 ATeV. The simulation is car- 
ried out using a microscopic simulator [6] incorporated 
in ALIROOT which is a GEANT3.21 and ROOT based 
simulation package used by the ALICE collaboration [7] . 

While we refer to [6] for deatil, in the following we 
briefly mention the salient features of the simulation. 
The ionization in the gas proceeds in two stages. Firstly, 
the electromagnetic interactions of the primary particles 
with the TPC gas (90% Ne and 10% C0 2 ) lead to the 
release of primary electrons with a statistics that follows 
a Poisson distribution. Thus, the distance S between two 
successive collisons leading to primary ionization can be 
simulated through an exponential distribution given by 
exp(-S/D)/D where D = {N prim f)- X is the mean dis- 
tance between primary ionizations, N pr i m is the number 
of primary electrons per cm produced by a Minimum Ion- 
izing Particle (MIP) and / is the Bethe-Bloch curve. At 
sufficient kinetic energy, the primary electrons produce 
secondaries creating an electron cluster with a total num- 
ber of electrons given by N to t — {Etot — Ipot)Wi + l where 
E to t is the energy loss in a given collision, W% is the effec- 
tive energy required to produce an electron-ion pair and 
Ipot is the first ionization potential. The electron clus- 
ter which is assumed to be point like undergoes diffusion 
while drifting towards the end-cap which is described by 
a three dimensional Gaussian with widths St — Dt\L~2 
and Sl = E>l\[Tm where Dt and Dt, are transverse and 
longitudinal diffusion constants and Ld is the total drift 
distance. An electron arriving at the anode wire creates 
an avalanche which induces a charge on the pad plane. 
The time signal is obtained by folding the avalanche with 
the shaping function of the pre-amplifier/shaper with a 
shaping time ~ 200 ns which is a compromise between 
the need for achieving a high signal to noise ratio and 
for avoiding overlap of successive signals. This signal is 
sampled with a frequency ~ 5.66 MHz which divides 
the total drift time of 88 fxs into about 500 time bins. 
The microscopic simulator also takes into account the 
electron loss in the drift gas due to presence of electron 
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negative gas like O2 and also ExB effect near the an- 
ode wires. The signal is digitized by a 10 bit A/D con- 
verter that generates Gaussian random noise with r.m.s. 
about 1000 e. Finally, the digitized data is processed and 
formatted by an Application Specific Integrated Circuit 
(ASIC) called ALTRO (ALICE TPC Read Out) [8]. A 
few typical parameters which are used in the simulations 
are given in table I. Although, we concentrate only on the 
TPC data, all component of the ALICE detectors as well 
as all passive materials are included in the simulation so 
as to create a situation close to the real experiment. 
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FIG. 1. (a) The amplitude versus time bin corresponding 
to a single pad. (b) An expanded view over 100 time bins. 

In ALTRO data format, zero suppressed data is 
recorded for each pad over all the time bins. This means, 
if we call bunch a group of adjacent over threshold sam- 
ples coming from one pad, the signal can be represented 
bunch by bunch. Figure 1 shows a typical plot of time bin 
versus amplitude for a given pad. Since the data is zero 
suppressed, it is sufficient to record three types of data, 
the sample amplitude in a given bunch, the bunch length 
and the time gap between two consecutive bunches. The 
amplitude distribution gives the energy loss spectrum 
with a long Landau tail and is not illuminating for the 
present purpose. Therefore, in the following, we will con- 
sider only the time gap and bunch length distributions 
built over all the pads. 

The time gap distribution corresponds to the distance 
between two trajectories along the drift direction (say 
z-direction). Figure 2 shows the plot of time gap dis- 
tribution at different multiplicities (M = 20000, 40000, 
60000 and 80000). Although, the tail of the distribu- 
tion is linear (in the semi-log scale) and depends on M 
as expected, it deviates from the linearity at the shorter 
distances. The above behavior can be described by a 
distribution of the type, 

f(z) = Az- f3 e- az for z>\ (1) 

The constant A is fixed by the requirement of normaliza- 



tion which gives A — [Lip(e a )] 1 where Li n (x) is the 
n th polylogarithm of x given by, 



Li n (x) = £ J 



(2) 



Defining z = Y^T=i z f( z )/J2T=i /( 2 )> the average of the 
above distribution is 



Li {e- a ) 



(3) 



Note that for — ► 0, the above distribution becomes a 
pure exponential with z = (1 — e _Q ) _1 . 
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FIG. 2. The frequency versus time gap for M = 20000, 
40000, 60000 and 80000 (top to bottom). The solid curves 
are the best fit obtained using Eq.(l) with different a and 
values as shown in figure 3. 

Figure 3 shows the parameters (filled circles in the 
left) and a (filled circles on the right) extracted from 
fitting the simulated data points. Note that a is a sta- 
tistical parameter which has a near linear dependence 
on the multiplicty M. However, the index increases 
with M and becomes a (nearly) constant at higer M val- 
ues. We will show below that this index is responsible 
for the dynamical scaling which does not depend on ex- 
ternal parameters like multiplicity M or magnetic field 
B above a critical value but depends on certain intrin- 
sic TPC parameters. Another interesting observation is 
that the product u — za also changes very slowly with M 
above a critical value of M <~ 30000. In order to appre- 
ciate the effect of dynamical scaling, we can plot figure 
2 in a reduced scale by dividing z by z and multiplying 
/ by z/A in Eq.(l). By this rescaling, the index will 
not change, but a will be rescaled to say u> = az. Since 
and oj are (approximately) independent of M above a 
certain value, the rescaled distributions will merge with 
each other as shown in figure 4. The solid curve is the fit 



of the type f(z) 



with ~ 0.45 and to ~ 0.67. 



Recall here that this scaling is quite similar to the famous 
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KNO scaling [9] which says that at high energies s, the 
probability distributions P n (s) of producing n particles 
in a certain collision process should exhibit the scaling re- 
lation nP n {s) = f(n/n). This scaling hypothesis asserts 
that if we rescale P n (s) measured at different energies via 
streching (shrinking) the vertical (horizontal) axes by n, 
the rescaled curves will coincide with each other. In the 
present case, the observed scaling is identical to KNO 
scaling if we replace energy s by the multiplicity M and 
n by average distance z although both corresponds to 
two different physical situations. 
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FIG. 3. The dynamical index /3 and u> = za (a) and a (b) 
as a function of Multiplicities in units of 1000. The solid line 
is the linaer fit a = 0.000000602M + 0.0104. 

Apart from this striking similarity with KNO, this 
rescaling is an elegant way to remove the statistical de- 
pendency from the dynamical behavior. It may be men- 
tioned here that such type of scaling is also expected in 
case of an exponential distribution when the exponent a 
is small . Under such limit (a « 1), since the average 
of an exponential distribution z = (1 — e _Q ) _1 <~ or 1 , 
lu = olz ~ 1 (The lu value may become more than unity if 
a is large). However, the interesting aspect of the present 
scaling is the deviation from unity and also having same 
u at all M. This means lu deviates from unity due to 
presence of the dynamical exponent (3. In absence of (3, 
the scaling would have followed the linear behaviour in 
the semi-log scale as shown by the dashed curve in fig- 
ure 4. Therefore, both (3 and r\ = 1 — u> arc dynamical 
exponents whose values do not depend on M. 

In the following, we investigate the origin of this dy- 
namical behavior and also the parameters that affect its 
values. The origin of this dynamical phenomena can 
be associated with the set of measurements which are 
strongly correlated. Although, at a given interaction 
point, the creation of primary electrons due to ioniza- 
tion is a random statistical process, the set of interac- 
tion points created by a single charged particles passing 
through the TPC are well correlated. Obviously, this cor- 



relation will depend on the multiplicty M and also on the 
number of measurements (number of pads, size etc). De- 
pending on these parameters, a critical value is reached 
beyond which the correlation becomes independent of M. 



10 r 




30 35 40 



FIG. 4. The frequency distribution in the reduced scale i.e. 
ujf/A versus z/z for different multiplicities. The solid line is 
a common fit with /3 = 0.45 and w = 0.67. The dashed line 
represents the exponential exp(— loz) with lu = 0.67. 

Although, the dynamical behavior may depend on 
other geometrical TPC parameters, in this study, we do 
not intend to change any such parameters as they have 
been optimized based on certain physics criteria. How- 
ever, we can see the effect of other relatively softer pa- 
rameters like diffusion constants, drift velocity, noise etc 
which are likely to change during operation. It is no- 
ticed that out of all dynamical parameters, r\ or (1 — ui) 
has some dependence on Dt and strong dependence on 
electronic noise. Since the produced electron clouds are 
broadened in transverse direction due to increased Dt, 
the correlation effect is also enhanced. It is found that 
increasing Dt by a factor of two, lu reduces only from 
0.67 to 0.62. Further increase in Dt has very little effect 
on lu. 

On the other hand, noise above the threshold reduces 
the correlation effect as it acts as an additional random 
source of electrons. As shown in figure 5, the middle 
curve is obtained with the parameters as given in the 
tablei (reffered as default parameters), while the upper 
curve corresponds to Dt twice the default value and the 
lower curve is due to an increased noise level from 1000 e 
to 2000 e. The correlation effect is lost due to noisy sig- 
nal for which lu — > 1. With increased noise level lu rises 
above unity (slowly) which is typical of an exponential 
behavior. Note that the dynamical parameters like Dt, 
Dt and drift velocity etc are expected to vary within 
< 4% during operation. This small fluctuation does not 
effect the lo value. On the otherhand, noise is quite un- 
predictable and also may go up due to the exposure of 
the electronics to the radiation for an extended period. 
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Therefore, the parameter u> can serve as an excellent on- 
line tool to monitor the noise level. 



fitting procedures. 
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FIG. 5. (a) Same as figure 4. The middle curve is with 
default parameters as given in the table. The upper curve 
is with Dt value increased by twice the default value where 
as the lower curve corresponds to increased noise level from 
lOOOe to 2000e for which w -> 1. 

So far, we have discussed only about the time gap dis- 
tribution. The tail of the bunch length distribution also 
shows a powerlaw behaviour ~ A n~ 7 as shown in fig- 
ure 6(a). Since the tail of the bunch length distribution 
corresponds to low energy electrons and delta rays, the 
exponent 7 is quite sensitive to the applied magnetic field 
(up to some limit). Figure 6(b) shows the bunch length 
distributions at different B values. May be this exponent 
can be used to monitor the magnetic field setting during 
the operation. Apart from this, we have not found the de- 
pendency of 7 on any other dynamical TPC parameters 
which could have been utilized for online monitoring. 

We would like to add here that the powerlaw exponent 
that affect the time gap distribution corresponds to a re- 
gion of short distances. Therefore, it is very difficult to 
extract the [3 parameters acurately and also fitting the 
time gap distribution with other functional forms like two 
exponentials can not be ruled out. However, we have seen 
that the quality of fit is much better with powerlaw ex- 
ponent. Further, as we have argued before, dynamical 
phenomena is expected in side a TPC due to many cor- 
related interaction points. The tail of the bunch length 
distribution having a perfect powerlaw behavior reflects 
this aspect rather unambiguously. Due to the same dy- 
namical origin, it is reasonable to assume a powerlaw 
with exponential cutoff for time gap distribution as well. 
We would also like to add here that the KNO type of 
scaling found in case of time gap distribution is indepen- 
dent of any fitting procedure and the deviation of the 
slope from unity is an indication that the noise level has 
remained below the threshold. This we consider as an 
important observation of this analyses irrespective of the 



FIG. 6. (a) The bunch length distribution in the reduced 
scale for different M values. The solid curve is a powerlaw fit 
with 7 ~ 4.8. (b) Same as (a), but for different B values. 

In conclusion, the time gap distribution shows an 
universal scaling behavior (KNO type) (at sufficiently 
large multiplicity) with an exponent deviating from unity. 
This is an interesting aspect to study the correlation phe- 
nomena in gas detector and also to understand physics 
(more explicitly, the physics that has gone into simula- 
tion) at the device level. An important pratical utility of 
this phenomena is the utilization of the above scaling ex- 
ponent to monitor the noise level above a given threshold 
without any computational complexity and also without 
building any rigorous models. In that sense, this analyses 
provide a model independent way to monitor the quality 
of the data what is being recorded. 



TABLE I. A few typical parameters used in the simulation. 
More details are given in ref [5]. 



Paramctrs 


Value 


Diffusion Constants {Dl — Dt) 


220 fim/^cm 


Drift Velocity at 400 V/cm 


2.83 cm/ (is 


Shaping Time 


190 ns FWHM 


Sampling Time 


200 ns 


Noise 


1000 c 


Magnetic field 


0.2 T 


Oxygen Content 


5 ppm 
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